Instability of black holes in massive gravity 
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We show that linear perturbations around the simplest black hole solution of massive bi-gravity 
theories, the bi-Schwarzschild solution, exhibit an unstable mode featuring the Gregory-Laflamme 
instability of higher-dimensional black strings. The result is obtained for the massive gravity theory 
which is free from the Boulware-Deser ghost, as well as for its extension with two dynamical metrics. 
These results may indicate that static black holes in massive gravity do not exist. For graviton mass 
of order of Hubble scale, however, the instability timescale is of order of the Hubble time. 



Massive gravity is an alternative to General Relativ- 
ity (GR) in which the graviton is supposed to gain mass. 
There are several motivations behind the study of mas- 
sive gravity: theoretical curiosity to construct a gravita- 
tion theory with a massive graviton, to add a new model 
in the zoo of alternative theories in order to test better 
GR compared to alternatives, explain present day cos- 
mological acceleration of the universe in a more "natu- 
ral" way than introduction of the cosmological constant, 
resolve the vacuum energy problem by hoping that the 
graviton becomes weaker a la Yukawa, etc. 

Massive gravity has started from the works by Fierz 
and Pauli [l| , who studied Lorentz invariant massive spin 
2 theory (hereafter referred to as FP), described by a 
quadratic action. In contrast to GR this model has five 
propagating degrees of freedom and it is excluded by So- 
lar system experiments due to an extra force associated 
with the helicity-0 mode. This extra force does not vanish 
in the limit of zero graviton mass, which is a consequence 
of the van Dam-Veltman-Zakharov discontinuity [2]. It 
was, however, conjectured by Vainshtein in [3| that GR 
can be restored at nonlinear level for nonlinear massive 
gravity theories. A nonlinear completion of FP can be 
achieved by introducing a second metric [4|, which is 
needed to construct a mass term (in the original linear FP 
model the role of the second metric is played by the flat 
metric coinciding with the background solution for the 
first metric). Whether the Vainshtein mechanism indeed 
works, namely whether there are solutions around mat- 
ter sources approximating GR inside some radius (the 
Vainshtein radius) that continuously join an asymptoti- 
cally flat solution of the linearized FP, was an unresolved 
question for a long time. Only recently it was shown that 
for a large class of nonlinear FP models the Vainshtein 
mechanism indeed works [5(. 

Another problem of massive gravity is the appearance 
of a ghost propagating degree of freedom at nonlinear 
level in generic nonlinear FP theory [6], the so called 
Boulware-Deser ghost. This problem has been also re- 
solved only very recently by de Rham, Gabadadze and 



Tolley 7] (dRGT theory), who found a two-parameter 
family of nonlinear FP theories free from the Boulware- 
Deser ghost. Similar to the general nonlinear FP 
model, the Vainshtein mechanism also operates in this 
model [a |Sj. The dRGT model can be generalized by 



promoting the second metric to be dynamical 12j. Al- 
though the dRGT theory and its extension are disputed 
from different angles 13l4l6j , these are the only candidate 
to be pathology-free nonlinear FP theory so far. There- 
fore, it is certainly worthwhile to study their different 
applications and physical consequences. 

In this paper we show that the simplest black hole 
solutions in both dRGT theory and its extension are 
unstable. Namely, the bi-gravity solutions with two 
Schwarzschild metrics simultaneously diagonal exhibit an 
unstable mode. The result is valid for bi-Schwarzschild 
solutions with both equal and proportional metrics. 

One can divide black hole solutions in bi-gravity theo- 
ries into two types. The first one corresponds to solutions 
in which the metrics are not simultaneously bi-diagonal, 
i.e. if one in some coordinates is diagonal, the other 
metric is not. It is possible to find some classes of solu- 
tions analytically, }9l4llj [see also recent reviews [l7j and 
references within]. Such solutions typically have a very 
non-trivial global structure, similar to those in the gen- 
eral nonlinear FP theory [18| and it is not clear if they are 
of physical relevance, in particular it may be difficult to 
form this type of solutions during gravitational collapse 
(if the result of gravitational collapse in massive gravity 
is a black hole). 

On the other hand, spherically symmetric static bi- 
diagonal solutions, which seem to be more natural out- 
come of gravitational collapse (at least from configura- 
tions with smooth sources featuring the Vainshtein be- 
havior), must share the same Killing horizon [19| |, if one 
assumes non-singular solutions. Therefore the variety 
of bi-diagonal solutions is highly restricted if the non- 
singularity of both metrics is assumed. 

It is, however, not difficult to extend the known so- 
lutions in GR to bi-gravity, by simply identifying the 



two metics with a solution for GR. This ensures that 
the interaction mass term is zero and the equations of 
motion are satisfied for both the dRGT model and its 
extension. For example, the analog of the simplest black 
hole solution in GR, the Schwarzschild solution, is the 
bi-Schwarzschild solution in bi-gravity. In this paper we 
will show explicitly that its linear perturbations satisfy 
the same equations studied by Gregory and Laflammc 
(GL) for the relevant four dimensional perturbations of 
higher-dimensional black string solutions |20|, featuring 
the well known GL instability. 

The action for the extension of the dRGT model with 
two dynamical metrics reads [7|, LL2| , 
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where Mp is the Planck mass, R g and IZf are the 
Einstcin-Hilbert terms for the g and / metrics corre- 
spondingly, k is a dimensionless parameter and U is the 
interaction mass term, 

U=U 2 +a 3 U 3 + a A U A . (2) 

The interaction can be written in terms of the matrix /C£J , 
K.% = 5% — \/g^ a jav 1 where the square root of a matrix 
A is understood as (VA)£(VA)" = A£. The pieces of 
the ghost-free interaction term <j2j) can be cast in the 
following form, 
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Varying the action with respect to g^ v and f^ gives cor- 
respondingly, in the vacuum, 



G^ = m 2 T" Q, 
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where Gp_ v and Q^ are the Einstein tensors for g^ v and 
ffj, v correspondingly, and Tj* v and Tff v come from the in- 
teraction term, 



Tjf v = -g,.? {€ - Ki) (K - K\S%) + o(/c 2 ), 
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where the higher order terms in K, we collectively denoted 
as C(/C 2 ). 

Note that the dRGT model for one dynamical metric 
is described by ^ without the last term, defining dy- 
namics of the metric /^„. In this case the metric f^ 
becomes non-dynamical and one should solve only the 
first equations in (|4]). 

If both metrics are identical, g^ v — / p „, then K, = 0, 
0. Therefore any vacuum GR 



and hence Tr v 



l±v 



solution g^ v is also a solution for the bi-gravity theory 
provided that /^ = g^ v . In particular, the simplest 
black hole solution in bi-gravity theory is obtained by 
identifying the line element for each of the metrics with 
the Schwarzschild one, 



ds l 
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where rg is the Schwarzschild radius. Now let us consider 
perturbations on top of the bi-Schwarzschild solution, so 
that we write 

g^v = g^ v + h^v, ffj, v — j ^ v + h^ v , (7) 

where g^J — fjiJ is given by ©. It is not difficult to see 
that g* a f av = 6»+ (K - h>t) and 
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where we introduced the notations h ' = M — h£ and 

u V V 

h = h£ + nh^. Linearization of (j4]) yields, 



£$h aP + ^(h$- gitv hM) =0, (9) 



£$h aP - ^ (h$ <?^ (_) ) = 
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where 
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+ D V - 9^0h + g^V a V f 3h a ^ + 2R\\h Xa ) 

is the linearized Einstein tensor on top of the metric g^ v 
in the vacuum, R^ v = R = 0. Combining equations (|9|) 
and (|10p we obtain, 

£$h$ + ^ (h& - ^ (_) ) = 0, (11) 
£$h$ = 0, (12) 

where m! = mWl + 1/k. Therefore hln, is a massive and 
hfu, is a massless spin-2 perturbation. Acting on (fTTj) 



with V" gives Vh^J = V M /i (_) and hence VV/iJ,; = 
\3h>~'. Using this last relation in the trace of (fTT|) one 
finds the conditions to be satisfied by h^: 

VftW = h^ = 0. (13) 

The equation for the massive spin-2 perturbation (|11[) 
with the conditions (1131) can be rewritten as 



Uh\z> + 2R° 



A h { - ] 



2 ht 



(14) 



Note that under the infinitesimal transformation associ- 
ated with the change of coordinates x M — > x M + £^ the 



massive spin-2 perturbation h\j, v does not change, while 
for the massless spin-2 perturbation h\u, this is the gauge 
transformation (which leaves the equations of motion un- 
changed) as for linearized GR. 

So far we considered the bi-dynamical extension of 
the dRGT model, where the two metrics are dynami- 
cal. In the context of the original dRGT model (with 
one dynamical metric, the other fixed and equal to 
Schwarzschild) , the same type of equation for the massive 
perturbations is valid. Indeed, in this case only perturba- 
tions of the dynamical metric are present, h^, and they 

obey Eq. (Ill|) with h]x U — > hp, v and n — > oo. 

The solutions to Equation (| 14[) with the conditions (|T5)) 
have already been studied in the context of the GL insta- 
bility 20]. In its simplest setting, linear perturbations 
around the five-dimensional black string 



(to a = -(l-^)* 9 +(l-^ 



dr 2 
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were considered. Performing a Fourier decomposition 
around the infinite 5th flat dimension, it was shown that 
the relevant perturbations are the 4D ones J dke lkz h\j,v 
(no instability was found in the h^z and h zz compo- 
nents), and that h^J satisfy the same massive spin two 
Equations (|T3J) with m' 2 = k 2 in the Schwarzschild back- 
ground. Unstable (Q > 0) spherically symmetric modes 
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regular at the future event horizon were found within the 
range 



< m! < 



O(l) 

r s 



(17) 



Based on these results we conclude that the bi- 
Schwarzschild massive gravity solutions are unstable too 
provided m' — my/l + 1/n satisfies (fTTj) . 

It is possible to extend these results to the case of pro- 
portional metrics, j^ v — cu 2 g^. By appropriate choice 
of «3, a4 and uj, one can obtain the bi-Schwarzschild so- 
lution for both the dRGT theory and its extension. The 
requirement of asymptotic flatness yields w as a function 
of «3 and «4 for the dRGT model, while for its exten- 
sion the asymptotic bi-flatness leaves only one-parameter 
family of solutions for to. One can check that the mass 
term in the perturbation equation(s) is merely modified 
by a factor depending on «3, a^ and co . The equations 



1 Note that in this case not only the linear terms in K in Eq. {5} but 
also the rest contribute to the mass term for the perturbations. 



for perturbations remain the same, up to the redefini- 
tion of the mass and therefore the instability is present 
in those solutions too if the inequality (ITT)) is satisfied 2 . 
Unfortunately, for non-bi-diagonal solutions it is not as 
straightforward to generalize our results. The main diffi- 
culty is that, unlike the case where the metrics are equal 
(or proportional), the matrix JC%, expressed in terms of 
h%, does not have a simple form like ((5$. Notice that 
f° r jp,v = ^ 2 giiu each element of K^ depends only on 
h£ (or h in the theory with two dynamical metrics) . 
This happens because in the case of proportional met- 
rics, computing Kjj boils down to taking a square root 
of (<S£ +h%), resulting in <5£ + \h^ + 0{h 2 ). Instead for 

the non-bi-diagonal case the expression 
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must be evaluated, where g^fa'J is not proportional to 
the unity matrix. An explicit analysis is needed to ad- 
dress the stability issue of such solutions. We leave the 
non-bi-diagonal case for future investigations. 

We would like to stress that the instability rate O = 
f2(m') depends on the mass of the graviton. As it was 
shown in [21|, the characteristic inverse time of the in- 
stability, for a fixed black hole mass M — M p rs/2, lin- 
early grows as a function of the graviton mass for small 



reaches a maximum at some m! 



and then de- 



m 

creases. In most situations that are physically relevant 
ra ~ m' « r s , since the graviton mass must be very 
small in order to satisfy solar system tests of gravity, 
say by the Vainshtein mechanism. In particular, to be 
responsible for the cosmic acceleration, one usually sets 
to ~ H, where H is the Hubble parameter. Therefore 
in the diagram of the instability [20] typically we are sit- 
ting in the regime of small graviton masses, which means 
that the characteristic instability scale is of order of the 
inverse graviton mass, r ~ mT 1 and is independent of 
the black hole mass. Therefore for 777 ~ H the character- 
istic time of the instability is of the order of Hubble time: 
for these scenarios, although the instability is present, it 
is not physically relevant for astrophysical black holes. 3 



It is interesting to discuss the recent work [22j in con- 
nection to our results, where black holes in the dRGT 
theory of massive gravity (with one dynamical metric, the 
other being fixed and equal to Minkowski) are considered. 
To get rid of the physical singularity of a bi-diagonal solu- 
tion, the authors assumed time dependence of the physi- 
cal metric, which results in the loss of mass by the black 
hole with the rate M ~ —r 2 s m 2 Mp. The "instability" 



2 There is a specific choice of parameters which gives zero mass 
term for the perturbation(s). This corresponds to a solution with 
infinitely strongly coupled perturbations. We do not consider 
this case here. 

3 We point out that in principle for small k (but bigger than 
(mrs) 2 to satisfy the inequality \VT)) the instability rate is en- 
hanced by the factor 1/ 



scale associated with the process is r ~ l/(rsm 2 ). This 
time scale is much larger than the characteristic time of 
instability we found, for rg <C to -1 . This means that the 
Gregory-Laflamme type of instability is much more effi- 
cient than the one found in 22[ . The difference in results 



must not be surprising, since both the physical processes 
leading to the instabilities and the black hole solutions 
studied in 22| and here are very different. The insta- 
bility of black holes found in [22j is due to the gradual 
flow of energy from the black hole, in fact it can be seen 
as an accretion of "phantom" energy onto the black hole 
with associated decreasing of the black hole mass [23j . 



The solution studied in |22j contains the flat reference 



metric, whereas we consider background solutions with 
two metrics being identical to the Schwarzschild metric. 
Moreover, the solution for the physical metric in [22| fea- 
tures the Yukawa decay far from the black hole. As it 
was pointed out in [22[ , the instability of such black holes 
was already conjectured in [24J. On the contrary, the bi- 
Schwarzschild solutions we study contain no asymptotic 
exponential decay and can be considered as solutions in 
the Vainshtein regime all the way from the horizon to 
infinity. 

There are several interesting questions that remain 
unanswered. First of all, the instability we discuss in this 
work operates only at linear level. It is not clear, how- 
ever, what would be the result of this instability after it 
grows to the nonlinear level, and, in general, what the 
fate of such unstable black holes is. It is worthwhile to 
mention that in the context of the black string instability 
in models with extra dimension(s), the linear instability 
grows to the nonlinear one, which results at the end in 
breaking up of the black string into an array of spherical 
black holes. The linear effective 4D description, however, 
ceases to work when the instability becomes nonlinear. 
Similarly, in bi-metric massive gravity we are not able to 
interpret the fate of bi-Schwarzschild black holes with- 
out the nonlinear analysis. It could be a stable black 
hole solution, if it exists, or somet hing else which is not 
a black hole. It was suggested in 22|, for a completely 
different physical process, that the end of that instability 
is Minkowski spacetime. On the other hand the insta- 
bility of the simplest (bi-Schwarzschild) black holes may 
indicate that static black holes in massive gravity do not 
exist and/or that they do not form in the process of grav- 
itational collapse. 

It would be also interesting to check if our results may 
be applied to other black hole solutions, in particular 
to non-bi-diagonal solutions, and in other viable mas- 
sive gravity theories, in particular in the Lorentz violat- 
ing massive gravity theory with 5 propagating degrees of 
freedom, recently proposed in 25 [. 

To summarize, we studied perturbations on top of 
black hole solutions in dRGT massive gravity theory and 
its bi-gravity extension with two dynamical metrics. We 
found that the simplest black hole solutions, namely the 



bi-Schwarzschild solutions, are unstable. Indeed, the per- 
turbation equations for such solutions coincide with those 
for the relevant perturbations of the higher-dimensional 
black string solutions featuring the well known Gregory- 
Laflamme instability. The instability time scale is of or- 
der of the graviton mass, therefore if the massive gravity 
is responsible for Dark Energy, then the instability is very 
mild, of the order of Hubble time. 
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Note added: After our submission, another paper [261 ] 
appeared in arXiv which confirms our conclusions on in- 
stability of black holes in massive gravity. 
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